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Abstract 

The impact of multiple ionization of air plasma on the 
inverse Bremsstrahlung absorption of laser radiation is 
investigated for air breathing laser propulsion. Thermo- 
chemical properties of multiply ionized air plasma 
species are computed for temperatures up to 200,000 deg 
K, using hydrogenic approximation of the electronic 
partition function; And those for neutral air molecules are 
also updated for temperatures up to 50,000 deg K, using 
available literature data. Three formulas for absorption 
are calculated and a general formula is recommended for 
multiple ionization absorption calculation. The plasma 
composition required for absorption calculation is 
obtained by increasing the degree of ionization 
sequentially, up to quadruple ionization, with a series of 
thermal equilibrium computations. The calculated 
second ionization absorption coefficient agrees 
reasonably well with that of available data. The 
importance of multiple ionization modeling is 
demonstrated with the finding that area under the 
quadruple ionization curve of absorption is found to be 
twice that of single ionization. The effort of this work is 
beneficial to the computational plasma aerodynamics 
modeling of laser lightcraft performance. 

Nomenclature 
a<> first Bohr radius 

aj - a 7 coefficients for thermodynamic functions 
bi, b 2 thermodynamic function integration constants 

c light velocity in vacuum 

C p molar heat capacity at constant pressure 
e proton charge 

E ionization potential 

Ei lowering of the ionization potential 


G Gaunt factor 

h Plank Constant 

I laser intensity 

H molar enthalpy at temperature for standard state 

Ho molar enthalpy at 0 K for standard state 

me electron mass 

n nth state of excitation energy 

n*,ni number density of electron and ions 

n a ,n ra number density of atoms and molecules 

Q partition function 

p pressure 

R universal gas constant 

S entropy at temperature for standard state 

s c Sackur-Tetrode constant 

si optical path length of the ith ray 

T temperature 

Z ion charge 


Creek 

Symbols 

a 

fine structure constant 

K 

absorption coefficient 

Kb 

Boltzmann’s constant 

E0 

permittivity of free space 

h 

Plank constant / 2 n 


real refractive index 

CO 

angular frequency 

Subscripts 

e 

electron 

i 

ith state 


Introduction 

Kantrowitz kl first suggested a new possibility for 
dramatic cost reductions in mass launching to Earth orbit 
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with a ground-based high-power laser in the 70’s. Since 
then, a propulsion system supported by a laser-sustained 
plasma has been the subject of many researches P1,G 13 1,62 
The main advantage gained by laser propulsion over 
chemical propulsion is the low- weight system obtained 
from decoupling the energy source from the vehicle, and 
high specific impulse resulting in low fuel consumption. 
In addition, the flame temperature of a combustion 
process is limited, whereas the propellant temperature 
reachable during laser propulsion can be several order- 
of-magnitude higher. 

Several air-breathing laser propulsion concepts have 
been demonstrated in the past few years. For example, 
recent publications show that a spin stabilized Myrabo 
lightcraft reached 71 m in record height during vertical 
free flights outdoors , 1 m while a different parabolic flyer 
(design) was propelled from the ground of the laboratory 
to its 8 m high ceiling. 62 High energy CO 2 lasers were 
used in both tests. 

Researches using computational plasma aerodynamics 
have also been making progress in the field of laser 
propulsion. For example, Molvik et al. M2 considered the 
interaction between a continuous laser beam and a 
flowing hydrogen gas using a structured- grid formulation 
and constant absorptivity. Jeng and Keefer J1 did similar 
analysis with an expression for the absorption coefficient 
at CO 2 laser wavelength of 10.6 pm considering both 
electron-ion and electron-neutral inverse Bremsstrahlung. 
Conrad, et al. CI modeled a continuous optical discharge 
stabilized by nitrogen gas flows in weakly focused laser 
beam, using a absorption coefficient formula at 10.6 pm 
ignoring second ionization of atoms. 01 Recently, Wang 
et al. performed transient performance calculations^ on 
a Myrabo lightcraft (energized by a pulsed laser beam) 
using an unstructured-grid formulation and the same 
single ionization absorption formula used by Conrad et 
al. 

In a computational plasma aerodynamics study, using 
the modeling of a Myrabo lightcraft as an example, the 
focusing of the laser radiation is solved first, followed by 
computing the initial air breakdown and the creation of 
seed free electrons. When enough seed electrons are 
produced, they absorb more photons and resulting in 
more air breakdown and producing more free electrons. 
An avalanche of free electrons soon follows and a strong 
shock wave is generated. These are all solved with 
transport equations of continuity, energy, momentum, 
and species continuity, along with physical models such 
as finite-rate chemistry, high temperature 
thermodynamics, beam attenuation through absorption, 
beam refraction, and non-equilibrium radiation. The 
computational model then computes the subsequent 
traveling of the shock wave through air, the heating and 
ionizing of which (such that the air plasma becomes 


capable of absorbing more laser radiation), and most 
importantly of all, the thrust of the lightcraft propelled by 
the shock wave. Figure 1 shows a snapshot of the 
computed (heavy gas) temperature contours and laser 
beam traces at an elapsed time of 18 ps. wl It can be seen 
that the laser beam reflects specularly on the optical 
surface and focuses onto a focal “point” on the shroud 
where the breakdown of air starts. The temperature 
contours in Fig. 1 also describes the growth of the plasma 
front. The “protrusion” of the plasma front indicates that 
the plasma front (and the shock wave) is propagating up 
the beam - a result of successive heating and ionizing of 
the medium (air) such that the medium becomes capable 
of absorbing more laser energy and propagates further. 
Figure 1 is an indication of the potential ability of the 
computational plasma aerodynamics in describing the 
optical breakdown phenomenon associated with laser 
propulsion. It also indicates the importance of a realistic 
absorption model since the propulsion physics start with 
the absorption of laser energy. It goes without saying the 
accuracy of the absorption model affects that of the 
computed performance. 

Other than the constant absorption coefficient used by 

Molvik et al., M2 all the variable absorption formulas used 

j 1 c 1 w i c 

by the afore-mentioned modeling efforts ’ ' are of 

the single ionization category. While applying single 
ionization formula to a hydrogen plasma J1 is reasonable, 
there is room for improvement when it is applied to 
nitrogen and air plasmas, 01 * Wl since the atomic numbers 
of nitrogen and oxygen atoms are eight and nine, 
respectively. Given that the single ionization formula 
was formulated by Raizer and Tybulewicz R1 in the 70 ’s, 
while a simplified procedure for calculating an averaged 
degree of multiple ionization was reported by Zefdovich 
and Raizer 21 in the 60*s, and the same single ionization 
formula was still used in the 90’s, ci led the author to 
speculate that the difficulty resides in the scarcity of 
reliable high temperature thermodynamic properties for 
the multiply ionized air atoms. To understand the impact 
of multiply ionized atoms on the absorption of laser 
radiation, reliable high temperature thermodynamic 
properties for the multiply ionized air atoms have to be 
developed, and that is the motivation of this study. 

Thermophysics C haracterization nf Multiply 
Ionized air plasma 

Thermo-Chemical System 

In Ref. Wl, the initial free electrons for plasma 
ignition and the subsequent avalanche of free electrons 
necessary for the optical breakdown were generated 
through the non-equilibrium, finite-rate air breakdown 
chemistry sub-model, where Park’s multitemperature air 
chemistry mechanism 63 was used. This mechanism 
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composes of the dissociation, NO exchange, associative 
ionization, charge exchange, electron impact ionization, 
and radiative recombination reactions. The eleven air 
plasma species used in this mechanism defines the 
thermo-chemical system for single ionization 
environment: N 2 , O 2 , NO, NO + , N, N + , O, 0 + , N 2 + , 02 + 
and e\ N 2 > O 2 , and NO are neutral molecules; N and O 
are neutral atoms; while NO + , N + , 0 + , N 2 + , 02 + are single 
ions. In order to consider multiply ionized air plasma 
atoms, up to quadruple ionization, six additional ions 
N* 2 , 0 +2 , N*\ O* 3 , N* 4 , and O* 4 must be added and their 
thermodynamic properties must be characterized. 

Hydrogenic Approxim ation of the Partition Function 
The high temperature thermodynamic properties of the 
six additional ions N +2 , 0 +2 , N^ 3 , O' 1 ’ 3 , In , and O* 4 can be 
expressed in terms of partition functions, following those 
formulated for monatomic gases. 02 For example, 

Cp 2 d 2 \nQ ^ dj^Q) ^5 
R dT 2 dT 2 

H_-Hp = J , d{ Ing) t 5 
RT dT 2 
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is chosen to be the energy of formation of the element 
from its reference state as defined in Gordon and 
McBride. 02 The result of the computation is a table of 
thermodynamic properties for different atoms and 
multiply ionized ions as a function of temperature and 
pressures. Pressure is calculated from the ideal gas 
equation of state plus the Coulomb pressure correction. 

Thermodynamic Function Generation 

The next step is to construct the three thermodynamic 
functions of heat capacity, enthalpy and entropy as 
functions of temperature in a form compatible to most 
computational plasma aerodynamics codes. The standard 
form of Gordon and McBride 02 is used: 


~~- = 0[T ^ + a 2 T ^ + <23 + a^T + a$T^ + + a~jT* 

JL = - axT “ 2 + a 2 T~ l lnT + a 3 + a 4 j + a 5 ~^ + a 6 ^ 


5 

-2 


+ an + - 

' 5 7 


S T~^ _i 7 2 7 3 

— = -ai — 1 +^3lnr + o 4 r + o 5 — +06 — 


+ 07 — + />2 


5 T d(\nQ) 
R dT 


+ ln<2 + — \nM + 


—In T + s r 

2 c 


In this study, the partition functions of the multiply 
ionized atoms are characterized with the method of 
hydrogenic approximation ZU2 That is, the multiply 
ionized atoms are treated as hydrogen-like atoms, 
represented by a system consisting of a positive nucleus 
with a charge Z and a single electron. The transformed 
electronic partition takes the expression 


Qi = £*2n 2 exp 






/J 


These coefficients are obtained through a least-square 
curve- fit procedure for each temperature interval. 
Following Gordon and McBride 02 , three temperature 
intervals are used in this study. Unlike Gordon and 
McBride 02 , the three temperature intervals are made 
different for different species in order to achieve best fit 
of the generated thermodynamic properties, since the 
heat capacities of the multiply ionized atoms peak at 
vastly different temperatures. To construct the enthalpy 
curve, the heat of formation of multiply ionized ions at 
reference state needs to be estimated. This is 
accomplished by writing an ionization reaction, e.g., for 
N + : 


where the summation is truncated when [Ei(l - 1/n 2 )] is 
greater than [Ej - ZEJ. The lowering of the ionization 
potential is proportional to the square root of the sum of 
the electron density plus the ion densities times their 
charge squared divided by temperature. For consistency 
purpose, the high temperature thermodynamic properties 
of N, N + , O, 0 + are also characterized as the hydrogen- 
like atoms. 

The computational procedure is set up such that the 
internal energy and density are used as input. The 
internal energy of a given species is a function of the 
temperature, partition function, energy of ionization, and 
reference point energy. This reference is arbitrary, but it 


N->N+ +e 

The heat of reaction of this ionization reaction is the 
ionization potential. The heat of formation of N + takes 
the form 

For validation purpose, the calculated heat of 
formations for singly ionized N + and 0 + are comparable 
to those published in Gurvich et al. 03 In addition, the 
calculated entropy of formations using hydrogenic 
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approximation of the partition function for N, O, N + and 
0 + are also comparable to those published in Gurvich, et 
al. as well. 

Figure 2 shows the computed heat capacities and curve 
fits for N, N\ N* 2 , N* 3 , and N 44 , while Fig. 3 shows those 
for O, 0 + , 0 +2 , 0 + \ and O 44 The peak heat capacity 
increases with temperature as the number of electrons 
stripped increases. The sensible heat capacities of the 
multiply ionized air plasma species cover a temperature 
range from 10,000 deg K to approximately 200,000 deg. 
K. The seven coefficients polynomials fit the computed 
heat capacities reasonably well. 

The thermodynamic functions of the rest of the air 
plasma species (N 2 , O 2 , NO, NO + , N 2 4 , 02 + and e‘) can be 
found from Gordon and McBride 04 where the calculated 
data from Gurvich et al. 03 were curve-fitted. However, 
the applicable temperature range for these species were 
only calculated up to 20,000 deg K, as shown in Fig’s 4- 
6. This temperature range appears to be too low, in light 
of the computed heavy gas and electron temperatures can 
go as high as 500,000 deg K during optical breakdown of 
air inside the focused region of a laser lightcraft. wi In 
addition, heat capacity data of many species do not level 
off to a value at higher temperatures, indicating possible 
overprediction of the heat capacity when extrapolated 
beyond 20,000 deg K. Nevertheless, there is no problem 
with electron since its heat capacity is a constant at any 
temperature. For species NO + , N 2 4 , and O 2 4 , the low 
applicable temperature range does not present a problem 
either, since only trace amount of these species are 
produced at conditions of interest. For molecular species 
N 2 , O 2 , and NO, Balakrishnan reported correlations for 
specific heats up to 50,000 deg K. B2 However, the 
formulas used by Balakrishnan were criticized as 
inadequate/ 3 Jaffe’s calculated heat capacities for N 2 , O 2 , 
and NO were based on summations over all vibration- 
rotation energy levels for all known bound electronic 
states, and a scheme for the partitioning of the internal 
energy into vibrational, rotational and electronic 
contributions was presented which consistently accounts 
for the nonseparable nature of the various energy modes. 
Jaffe’s work appears creditable and is used in this study. 

Figures 4-6 show the heat capacities from all three 
sources for N 2 , O 2 , and NO, respectively. It can be seen 
that the heat capacities of Jaffe agree with those of 
Gordon and McBride reasonably well, while the heat 
capacities of Balakrishnan agree with those of the other 
two sources only at lower temperatures. The heat 
capacities of Jaffe are chosen for curve fitting for species 
N 2 and O 2 . For species NO, Gordon and McBride’s data 
were used for curve fitting up to 20,000 deg K, then 
Jaffe’s heat capacities were fitted for the higher 
temperatures. Note that all three curves approach their 
asymptote values near 50,000 deg K, meaning these 
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curves can be extrapolated to much higher temperatures, 
say, 500,000 deg K. It is a moot point though since these 
species do not survive beyond 20,000 deg K. 

leaser Absorption 

** -yy j 

In computational plasma aerodynamics modeling of 
the laser lightcraft flowfield where geometric optics is 
used to simulate the local intensity of the laser beam, the 
laser beam can be split into a number of individual rays. 
In the presence of absorption, the local intensity of each 
ray follows the Beer’s law: 


Through inverse Bremsstrahlung absorption, or free- 
free absorption, the rays are attenuated by free-electrons 
in its path. The three types of inverse Bremsstrahlung 
absorption, depending upon what kind of particle the 
electron is near when a photon is absorbed, are electron- 
ion, electron-atom, and electron molecular absorptions. 
According to Raizer and Tybulewicz, R1 the long- 
wavelength infrared radiation of a CO 2 laser at 10.6 pm 
is absorbed mainly by free electrons when it collide with 
ions. Hughes H1 gives a theoretical derivation of the 
electron-ion inverse Bremsstrahlung absorption 
coefficient for radiation at frequency to: 

w g n,-Z 2 e 6 g[l-exp(-ft<o/K: fi r)Y m g V ' 2 

03 fife ^chco 3 m]; \ (mk b T t 

The advantage of Hughes’ formula is its flexibility. 
That is, it can be used for radiation in any wavelength. In 
contrast, started with a different formulation from that of 
Hughes, corrected for stimulated emission in the single 
ionization range, assumed hay*CBT « 1, omitted the 
factor affecting only the photoionization, substituted hco 
for CO 2 laser wavelength, a formula of electron-ion 
inverse Bremsstrahlung absorption coefficient is 
expressed by Raizer and Tybulewicz: RI 


*C02 = 


10-4 pSg 

(7/ 10 4 ) 7/2 


where 


g = 0.551n 


27(7/10 4 )* /3 j 


- 1/3 


Note the above absorption coefficient formula does not 
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take the second ionization (or higher) into account. Also, 
electron pressure is used in lieu of electron number 
density. Both formulas described above consider only 
the electron-ion inverse Bremsstrahlung absorption. On 
the other hand, Mertogul M3 computed the absorption 
coefficient using all three types of inverse 
Bremsstrahlung. The formula used for the calculation of 
electron-ion inverse Bremsstrahlung absorption 

coefficient is that given by Stallcop: 51 



in which the free-free Gaunt factor was curve fitted form 
reported data of Karzas and Latter at wavelength of 10.6 
pm and T < 10,000 deg K: 

g = 1.07 + 6.9643jd0~ 5 r - 2.6786x10 ~V 2 

and for T > 10,000 deg K: 

g = 1.50 + 1. Oxio _5 r 

The expression used for the electron-atom inverse 
Bremsstrahlung absorption coefficient in the infrared 
limit is that given by Stallcop: 52 

K w =n e n a K B T^A5x\0~ 29 )^ 

exp(-4.862*7-(l- 0.2096/1:7- + 0.017*^ -0.00968*:^) 
k T 

and 



The expression used for the electron-molecule inverse 
Bremsstrahlung absorption coefficient is that given by 
Caledonia et al. 02 from the work of Dalgarno and Lane: 01 

*w =«,«„, (4.5 brlO -44 ) 

j-3/2 


where D is a complicated power series represented as a 
function of haVicsT and was given in a Appendix of 
Mertogul. M3 For a CO 2 laser wavelength of 10.6 pm, this 
expression of electron-molecule inverse Bremsstrahlung 
absorption coefficient is only valid for temperatures less 
than 4321.5 deg K. In addition, MertoguFs formula was 
derived for hydrogen plasma, a deviation from our 
interest in air breathing laser propulsion. Nevertheless, 
MertoguFs formula is included in this study to compare 
the relative importance of these three types of inverse 
Bremsstrahlung absorption. 

Results and Discussion 

Thermodynamic functions generated in this study for 
multiply ionized ions, atoms, and molecules are used as 
data base for a series of constant pressure ( 1 atmosphere) 
and temperature thermal equilibrium computations, in 
order to obtain the necessary compositions of electron, 
ions, atoms and molecules for laser absorption coefficient 
calculations. In the temperatures of interest, equilibrium 
state is probably a reasonable assumption. Minimization 
of free energy of a thermo-chemical system, similar to that 
described in Gordon and McBride, 04 is used as the 
algorithm for achieving the equilibrium state and is not 
repeated in here. Figure 7 shows the air plasma species 
compositions considering single ionization only. As 
temperature increases, the molecules disappear quickly and 
atoms emerge. And then atoms disappear, while electron 
and ions (N* and O*) rise, eventually, the species 
concentrations of electron and ions level off at about 
32,000 deg K. The final electron mole fraction of 0.5 is 
the result of single ionization. Note that the 
concentrations of NO + , N 2 + , and C> 2 + are indeed 
negligible. 

Figure 8 shows the comparison of calculated 
absorption coefficients of air for CO 2 laser radiation 
using the information from Fig. 7. Also plotted in Fig. 8 
are points read off from Fig. 6.18 of Raizer and 
Tybulewicz R1 while allowing for double ionization. At 
low temperatures, the calculated absorption coefficients are 
extremely low and rise sharply with increasing temperate 
around 8000 deg K. The rise then slows down and the 
absorption passes through a maximum. The calculated 
absorption drops monotonically from the maximum as the 
system completes the single ionization, eventually nearing 
zero absorption around 80,000 deg K. 

It can be seen that the curves using Raizer and 
Tybulewicz formula and that of Hughes are reasonably 
close, although the peak value of Hughes is higher. The 
peak value of Mertogul is the lowest among the three. 
This is not surprising since MertoguFs formula was 
meant for hydrogen plasma. Of interest is at lower 
temperatures where the electron-atom, and electron- 


1-exd 


hco ^ 


*bT\ 
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molecular inverse Bremsstrahlung absorptions should 
show some contribution, but none was observed, 
indicating the free-free inverse Bremsstrahlung 
absorption is indeed the main absorption process among 
the three as described by Raizer and Tybulewicz. 

On the other hand, the curve of Hughes Formula 
agrees better with data points from Fig. 6.18 of Raizer 
and Tybulewicz than that calculated using Raizer and 
Tybulewicz formula, until the second ionization takes 
place. This is somewhat perplexing and it is speculated 
that difference in generating the electron pressures may 
have been the problem. Raizer and Tybulewicz did not 
elaborate how the electron pressure was obtained, nor did 
they explain how the second ionization portion was 
arrived and only the single ionization formula was given. 
Nevertheless, the mathematical reason for the monotonic 
decrease of the calculated absorption coefficient after the 
passing of the maximum can be seen clearly from Fig. 7 
that electron partial pressure is a constant as temperature 
exceeds 30,000 deg K, that the Gaunt factor is essentially 
a constant since it ranges from 2.3 to 3.2 in the 
temperature range of interest, and that the denominator of 
Raizer and Tybulewicz formula eventually grows as 
temperature increases. The downward trend of all three 
curves after their peaks confirms the single ionization 
system does not produce a second rise in absorption 
coefficient. 

Although only reported to 27,000 deg K, Rl the 
importance of second ionization is evident from Fig. 8 in 
which the single ionization formulas under-predict the 
absorption coefficient at high temperatures where double 
ionization occurs. This has strong implications for laser 
lightcraft performance computations using computational 
plasma aerodynamics. The implication of double 
ionization also implies the potential importance of triple 
ionization, ..., etc. That raises an issue of rising 
computational cost if too high a degree of ionization is 
considered however, for the computational cost is 
proportional to the square of the number of species 
considered. 

Based on the result in Fig. 8, the general Hughes 
formula is used to investigate the effect of multiple 
ionization, in conjunction with the Gaunt factor given by 
Raizer and Tybulewicz. This is accomplished by 
performing a series of equilibrium computations using 
the characterized thermodynamic properties for multiply 
ionized air ions N 42 , 0 +2 , N +3 , O 43 , N 44 , and 0 44 . The 
plasma composition required for absorption calculation is 
obtained by increasing the degree of ionization 
sequentially, up to quadruple ionization. 

Figures 9-11 show the equilibrium air plasma species 
compositions for double, triple, and quadruple 
ionizations, respectively. As expected, the surviving ions 
for double ionization are N 42 and 0 42 for temperatures 
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exceeding 50,000 deg K, those for triple ionization are 
N 43 and O 43 for temperatures exceeding around 95,000 
deg K, and those for quadruple ionization are N 44 and 
O 44 for temperatures exceeding 100,000 deg K. The 
peaks of those multiply ionized species concentration 
decrease as the degree of ionization increases. Most 
importantly, the mole fractions of the free electron 
increase from 0.5 for single ionization, to 0.67, 0.75, and 
0.8 for double, triple, and quadruple ionizations, 
respectively. 

Figure 12 shows the comparison of electron number 
densities for double, triple, and quadruple ionizations. A 
second peak of the electron number density occurs at 
around 30,000 deg K, due to the double ionization. Triple 
ionization increases the overall electron number density 
from approximately 40,000 degree up. Quadruple 
ionization increases the electron number density from 
around 65,000 degree up, but the amount of increase 
becomes less as the degree of ionization increases. That 
indicates for computational purpose, quadruple ionization 
is probably enough. 

Figure 13 shows the comparison of calculated 
absorption coefficient curves for double, triple, and 
quadruple ionizations. A second maximum in the 
absorption coefficient at temperatures near 30,000 deg K 
occurs due to double ionization. As the system completes 
the second ionization, the absorption again passes a 
maximum, although it is much less obvious as the last 
peak, and so on. The earlier part of the second ionization 
curve agrees well with the available points from Raizer 
and Tybulewicz. Note that the area under the quadruple 
ionization curve is about twice that of the single 
ionization. In summary, under a single ionization 
system, the absorption of photons drops sharply beyond 

25.000 deg K and appears to cease absorbing energy at 

80.000 deg K. When allowing for quadruple ionization, 
the absorption not only passes through multiple 
maximums, but also continues to absorb energy beyond 

100.000 deg K. 

Conclusions 

A thermophysics characterization of inverse 
Bremsstrahlung absorption of laser radiation is performed. 
Thermo-chemical properties of multiple ionized air plasma 
species are generated using hydrogenic approximation of 
the electronic partition function and those for neutral air 
molecules are also generated using updated literature data. 
Three formulas for absorption are calculated and a general 
formula is recommended for multiple ionization absorption 
calculation. A series of thermal equilibrium computations 
are performed to show the effect of multiple ionization on 
the free electron concentration and on the inverse 
Bremsstrahlung absorption coefficient. The calculated 
second ionization absorption coefficient agrees 
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reasonably well with available data of literature. In 
addition, it is found that the area under the quadruple 
ionization curve of absorption is about twice that of the 
single ionization. The result of this study can be applied 
to the computational plasma aerodynamics modeling of 
laser propulsion physics. 
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T r deg K 

Fig. 2 Computed heat capacities and curve fits for 
multiply ionized nitrogen atoms N, N*, N 4-2 , N + \ and N 44 . 



temperature contours and laser ray traces for a Myrabo T K 

lightcraft at 18 pis. Contours scale: 0 - 24180. 

Fig. 3 Computed heat capacities and curve fits for 
multiply ionized oxygen atoms O, 0 + , 0 +2 , 0 + \ and 0 44 . 
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T, deg K 

Fig. 4 A comparison of the heat capacities and curve fit for 
neutral N 2 . 



T, deg K 


Fig. 5 A comparison of the heat capacities and curve fit for 
neutral molecule O 2 . 



T. deg K 


Fig. 6 A comparison of the heat capacities and curve fit for 
neutral molecule NO. 
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Fig. 7 A comparison of the equilibrium air plasma species 
compositions for single ionization. 
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Fig. 8 A comparison of the calculated absorption 
coefficients of air. 



T, deg K 

Fig. 9 A comparison of the equilibrium air plasma 
species compositions allowing for double ionization. 



T. deg K 

Fig. 10 A comparison of the equilibrium air plasma 
species compositions allowing for triple ionization. 



T. deg K 

Fig. 1 1 A comparison of the equilibrium air plasma 
species compositions allowing for quadruple ionizations. 
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Fig. 12 A comparison of the electron number densities 
for double, triple, and quadruple ionizations. 



Fig. 13 A comparison of calculated absorption coefficients 
for double, triple, and quadruple ionizations. 
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